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ABSTRACT 

We present distributions of colors of stars along the horizontal branch of the 
globular cluster NGC 2808, from Hubble Space Telescope WFPC2 imaging in B, V, 
and an ultraviolet filter (F218W). This cluster's HB is already known to be strongly 
bimodal, with approximately equal-sized HB populations widely separated in the 
color-magnitude diagram. Our images reveal a long blue tail with two gaps, for a 
total of four nearly distinct HB groups. These gaps are very narrow, corresponding 
to envelope-mass differences of only ~ O.OIMq. This remarkable multimodality may 
be a signature of mass-loss processes, subtle composition variations, or dynamical 
effects; we briefly summarize the possibilities. The existence of narrow gaps between 
distinct clumps on the HB presents a challenge for models that attempt to explain HB 
bimodality or other peculiar HB structures. 



Subject headings: globular clusters: individual (NGC 2808) — stars: horizontal branch 



1. Introduction 

Despite decades of work, the causes of the distribution of colors of horizontal-branch stars 
in globular clusters remain a mystery. In general, more-metal-rich clusters have red HBs, while 
metal-poor HBs are bluer. But there is a large degree of variation in HB morphologies, and a 
number of cluster HBs do not have the color expected for their metal abundance. The distribution 
of observed colors requires that up to ~ 0.3Mq have been lost by each star (Castellani & Renzini 
1968, Iben & Rood 1970). The mechanism for this mass loss is not understood, though, and our 
lack of understanding of this phenomenon pervades all discussion of this subject. 

A particularly strange piece of this puzzle has been the few clusters that have bimodal HB 
distributions, including NGC 1851 (Stetson 1981), 2808 (Harris 1974), and 6229 (Borissova et al. 
1997) at moderate metallicities. Perhaps even more surprising are the newly discovered bimodal 
HBs in the metal-rich clusters NGC 6388 and 6441 (Rich et al. 1997). A set of extremely hot HB 
stars is also found in the metal-rich open cluster NGC 6791 (Kaluzny & Udalski 1992, Liebert, 
Saffer, & Green 1994), whose HB is predominantly red. As well as these red-blue bimodalities, a 
number of clusters appear to have gaps along the blue HB tail, with the bluest stars distinguished 
from the slightly cooler stars by an obvious break in the HB sequence. 

In this paper we focus mainly on the blue HB of NGC 2808. We present new observations of 
the central part of the cluster, taken with the Wide Field/Planetary Camera 2 (WFPC2) aboard 
the Hubble Space Telescope (HST). These images contain a larger sample of post-main-sequence 
stars in the cluster than any taken previously. They show that the blue part of the cluster's 
HB is considerably more extended than previously known, and reveal two narrow gaps along the 
blue HB. The overall distribution of HB stars is thus multi-modal, with at least four significant 
populations. 

The distinct clumps along the blue HB tail, and the narrow gaps that separate them, present 
a special challenge for models of bimodal and other peculiar HB morphologies. 



2. Observations and Reduction 

NGC 2808 was observed by the WFPC2 on 2 December 1995. Several exposures were taken 
through the F555W and F439W filters, which approximate Johnson V and B; total exposure 
times were 57 s and 510 s, respectively. A longer series of exposures (totaling 3300 s) was taken 
through F218W, a UV filter that peaks near 2180A. The images were processed by the standard 
HST pipeline; we then used DAGPHOT (Stetson 1987) to identify stars and measure magnitudes. 

The photometry was calibrated to Johnson V and B (for F555W and F439W), and to the 
STMAG system (for F218W), using zeropoints and transformations given by Holtzman et al. 
(1995). The effects of contamination in the UV image were removed according to the prescription 
of Whitmore, Heyer, & Baggett (1996). Our calibration is ~ 0.1 mag fainter in V, and ^^ 0.1 mag 



bluer in B — V, than that of Ferraro et al. (1990); Ferraro (1996, private communication) has 
confirmed that his recahbration of their data is consistent with these differences. 

The {y,B — V) color-magnitude diagram (CMD) is shown in Figure 1, which contains over 
27,000 objects. The previously known portions of the HB are visible near (V = 16.3, B — V = 0.8) 
and (y = 17, B — V = Q.l). What is new in Fig. 1 is the extension of the blue tail of the HB, all 
the way down to F ~ 21. Moreover, even from this diagram it is clear that the blue tail is not 
uniformly populated: there are discontinuities in the density of points near y = 18 and V = 20. 

The structure of the blue HB becomes more apparent in the (i?,F218W — B) diagram in 
Figure 2, which contains 437 objects. The stars farthest to the left in the lower panel are at the 
bottom of the blue tail in Fig. 1. The stars near [B = 18, F218W — B = 1) correspond to the 
brightest blue stragglers in Fig. 1. The red HB clump and the RGB are faint in F218W, and 
do not appear in Fig. 2. The upper panel of Fig. 2 shows a histogram of F218W — B colors of 
HB stars. There are three peaks, at F218W — i? ~ —0.5, —1.7, and —2.4, and two gaps, near 
F218W — B = —1.9 and —1.3 (not entirely empty). When the red HB clump is included, the HB 
of NGC 2808 is now seen to be four groups of stars along a single sequence. 

We ran completeness experiments using the F218W and B images. The counts are virtually 
100% complete for HB stars brighter than B = 20, and are ~ 90% complete for the bluest clump, 
averaged over the observed area. The small amount of incompleteness will not dramatically change 
the relative heights of the three peaks of the histogram in Fig. 2, and certainly could not be the 
cause of the gaps. 

The four groups consist of ~ 350, 275, 70, and 60 stars, respectively, going from red to blue. 
The Lee index [B — R)/[B + V + R) is thus ^ 0.1, in contrast to previous studies that have found 
more red stars than blue (Byun & Lee 1993). 



3. Physical Parameters along the HB 

We fitted the zero-age HB (ZAHB) and evolutionary models of Dorman, Rood, & O'Connell 
(1993) to both observed CMDs. Newly computed models extend the ZAHB to envelope masses 
-^cnv = 10~ Mq. Stars with such a small envelope are essentially on the He main sequence, since 
H-burning is negligible throughout their HB evolution. Models computed with more modern 
physical inputs make no difference for this paper (cf. Yi, Lee & Demarque 1993, Hill et al. 
1996). We used Kurucz (1992) synthetic stellar fluxes to generate model colors, and assumed 
Ay21sw / E{B — V) = 8.74, using the extinction curve of Cardelli et al. (1989). The extension of 
the ZAHB redward of the red HB clump is an artifact of the Kurucz-model colors; colors produced 
by Bell (1997, private communication) do match the clump. 

We could not flnd a combination of (m — M)q and E{B — V) that flt both CMDs with the 
nominal filter zeropoints, probably due to uncertainties in the calibration and in the transformation 



from the theoretical to the observational plane. However, by shifting the zeropoints (within the 
uncertainties given by Holtzman et al), we find excellent fits with parameters ranging from 
[{m — M)q,E{B — V)] = [15.25,0.16] to [15.4,0.09], depending on how we divide the zeropoint 
shift between the optical color and the F218W magnitude. The canonical values are 14.8 and 0.22, 
respectively (Ferraro et al. 1990, Peterson 1993), but the new UV colors and revised V magnitudes 
rule out those values, even with the maximum possible shifts. Note that most of the difference 
between the old and new values is along the reddening vector. 

In Figure 3 we show the model fit to the data for one fit (15.25 and 0.16), with the HST B — V 
colors shifted 0.07 mag to the red. The HB evolutionary tracks shown are for models with masses 
of 0.540, 0.530, 0.510, 0.495, and O.4848M0, going from red to blue, with Mcorc = O.4847M0. The 
bluest track is for a star with virtually no envelope, and thus represents a theoretical limit to the 
HB. The isochrone shown is for a population with [Fe/H] = —1.48 and an age of 14 Gyr (not fitted 
to the data near the turnoff). The zero- age main sequence shown extends up to 1.5M0. 

The two blue HB gaps are at Tgg near 17, 000 and 25, 000 K, and masses near 0.54 and 
O.495M0, respectively. In both cases the gap widths correspond to mass-loss differences of only 
~ O.O1M0. (The other gap, between the red and blue HBs, is wider: ~ O.O5M0.) The two 
intermediate groups of stars are ^ O.O5M0 wide, but the bluest clump spans only O.OIM0 in mass 
(but nearly 10, OOOK in Tcfj). Whatever physical mechanism is responsible for the presence of 
the blue tail must also be capable of producing such astonishingly narrow gaps between discrete 
clumps. 

The uppermost gap on the blue HB is close to the mass where the HB evolutionary tracks 
change from predominantly redward to blueward, as a result of the decreasing contribution of the 
hydrogen-burning shell. This change may account for part of that gap, although it seems likely 
that a gap in the mass distribution exists as well. 



4. Discussion 

Why do narrow gaps appear on the horizontal branch of NGC 2808? Canonical stellar models 
cannot account for the gaps without appeal to the mass distribution (Whitney et al. 1997, and 
Fig. 3), with the possible exception of the middle gap. We must therefore consider the factors that 
determine a star's position on the ZAHB. 

It is not clear whether the red-blue bimodality necessarily has the same cause as the gaps 
within the blue HB. NGC 1851 shows the red-blue gap, but not the blue tail. On the other hand, 
there is significant variation in the morphologies of blue HBs. The HBs of NGC 6752 and M13 
(Paltrinieri et al. 1997) and M15 (Buonanno et al. 1985) show gaps that are arguably similar to 
the blue HB of NGC 2808. Two gaps have also been found in the blue HB of the field halo stars 
by Newell (1973). Thus, the bizarre appearance of the NGC 2808 HB may be the result of a 
combination of effects, some seen elsewhere in lesser degree, and some unique to this cluster. 



Variations in the overall HB morphology between clusters have received much attention. In 
addition to the "first parameter" (metal abundance), a number of "second-parameter" candidates 
have been proposed. Much recent debate has centered on whether the predominant second 
parameter is age (Lee, Demarque, & Zinn 1994; Stetson, VandenBerg, & Bolte 1996). Other 
candidates include He and CNO abundance (which alter the shell-burning rate), rotation and 
number density of stars (which, indirectly, affect the amount of mass lost by a star), or a 
combination of candidates. 

Most of these candidates do not offer a satisfactory explanation for the gaps we see here, or for 
red-blue bimodality in general. If NGC 2808 had a spread in [Fe/H] or age, we would not expect to 
see as narrow an RGB and turnoff, respectively, as are seen in the optical CMD (Buonanno et al. 
1984). The suggestion of van den Bergh (1996) that cluster mergers are responsible for bimodality 
appears unlikely for this same reason. A peculiar distribution of abundances of other elements 
could help to account for the HB structure seen here. However, CNO abundance variations have 
not been very successful in explaining the HB morphologies of other clusters (Dickens et al. 1991), 
and the problem of why the distribution was abnormal in NGC 2808 would remain. 

Various authors (Renzini 1977, Buonanno et al. 1985) have suggested that stellar rotation 
might account for bimodality. Rood et al. (1993) proposed that if rotation causes mass loss 
through some MHD effect, then there could be a critical rotational speed beyond which mass loss 
is enhanced. Another possibility is suggested by the hypothesis of Sweigart (1997) of He mixing 
on the RGB. Mixing would lead to bluer HB colors, because of the higher RGB tip luminosity and 
thus greater mass loss, and also because of the higher He abundance in the envelopes of the HB 
stars. If mixing took place only above a critical rotational speed, it might lead to a peculiar HB 
morphology. However, rapid rotation would lead to a larger core for the bluer stars, destroying 
the agreement of the HB locus and the constant-core-mass model prediction (Fig. 3). 

The local environment may play a role. Clusters with denser cores often have bluer HBs (Fusi 
Pecci et al. 1993, and references therein), and tidal stripping of red-giant envelopes during stellar 
encounters has been suggested as the cause. The revised parameters given in Fig. 3 lead to a lower 
central density for NGC 2808 than previously thought (logpo = 4.45, rather than 4.63 [Djorgovski 
1993], in Lqv/pc^). This density is still high, and the presence of the blue tail is consistent 
with what is seen in other dense clusters (Buonanno et al. 1997). We estimate that ~ 100 close 
encounters have occurred in the cluster in the last Gyr, although not all of these involved red 
giants. If tidal interactions lead to mass loss, then some stars would, by chance, have multiple 
encounters, which would explain the multiple HB populations in a natural way. 

One potential problem with this scenario is that the blue and red HB stars have identical 
radial distributions. Presumably, the progenitors of the blue HB stars suffered their encounters 
within the cluster core. A Fokker-Planck simulation of NGC 2808 shows that the orbits of stars 
initially within Tc do not evolve much in 10^ yr (the HB lifetime). The lack of radial gradients 
would thus require the encounters to occur near the base of the RGB. (The longer time spent 



at that stage would compensate for the stars' smaUer cross-section.) It is also possible that 
encounters involving binaries would expel stripped stars to larger radii. 

Also, for narrow gaps to appear on the HB, the encounters must be "tuned" to cause 
particular amounts of mass loss on each encounter — amounts that are nearly the same from star 
to star. It seems implausible that the sequence of tidal events would be so reproducible. Another 
objection is that the relative numbers of stars in the groups are not well fit by the predictions of 
this scenario, in its simplest form. If we assume that each star has an equal encounter probability 
per unit time (a poor approximation, given that the density varies widely with r), then the counts 
should follow a Poisson distribution. The best-fitting distribution predicts 352, 289, 119, and 40 
stars in the four groups, from red to blue. The last two predictions disagree strongly with the 
data: the x^ distribution gives a probability less than 10~^ of such a severe disagreement. On the 
other hand, a model that took the varying density into account might fit more closely. 

Finally, stellar encounters could not be the cause of the gaps in the field-star distribution, 
although the latter might arise from an entirely different phenomenon. The encounter frequency 
is also much lower in NGC 6229 and in M13. 

The presence of a binary companion could be important. A somewhat higher frequency of 
hot field subdwarfs are found in binaries (Allard et al. 1994), as are a number of the blue stars in 
NGC 6791 (Green et al. 1997). The bluest clump of HB stars in F218W - 5 has a shght redward 
extension in B — V ^ possibly indicating that some of these stars are blended with a main-sequence 
companion. However, the binary-star hypothesis would also lead us to expect radial population 
gradients, which are not observed. 

Since the physics of stellar mass loss is so poorly understood, we must be somewhat speculative 
in considering further options. One plausible hypothesis to explain multimodality is that multiple 
mass-loss processes act in red-giant stars. There is little physical justification available for any 
specific mass distribution, and no survey has addressed this problem in a systematic way. Some 
red HBs could even be multimodal, but this effect could not be observed because of the small 
variation of color with envelope mass for red HB stars. 

Mass-loss processes typically considered for solitary cool stars include molecular-line- and 
grain-driven winds, by which the material in the lower atmosphere must couple to the radiation 
field. Such processes might depend on the detailed line formation; if so, the locations of HB gaps 
may be consistent between clusters. Also, spectroscopy might reveal differences between stars on 
either side of a gap (see Peterson, Rood, & Crocker 1995, Moehler et al. 1995). 

In summary, we can provide no satisfactory explanation for the multimodality of the 
horizontal branch of NGC 2808. These observations serve to remind us that stellar mass loss is 
a poorly understood phenomenon that awaits much further work. Our hope is that this work 
will eventually tell us why NGC 2808, out of all well-observed globular clusters, is blessed with a 
combination of factors that gives its HB such a spectacular appearance. 
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Fig. 1.— (y,B - V) CMD of 27,286 stars in NGC 2808. Note the extended blue HB, and the 
prominent blue straggler sequence. 
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Fig. 2.— Lower panel: {B, F218W - B) CMD of 437 blue stars in NGC 2808. B magnitudes are in 
the Johnson system, and F218W magnitudes are in the STMAG system. The reddest portion of the 
HB does not appear on this diagram. The larger symbols are the HB stars that are histogrammed 
in the upper panel; smaller symbols are objects not on the HB. Upper panel: Histogram of colors 
of blue HB stars. Note the two gaps. 



13 



NGC S80e: (in-M^ = 15. E5, E(B-V) = 0,16 



- 



o 

as rn 

PI s 

O 



^4- 




0,5 

(B - V)^, [Johnson] 




-3 -3 -1 

(F218¥ [STMAG] - B [Johnson])(j 



Fig. 3. — A fit of ZAHB models and isochrones to the CMDs. Tlie upper panel shows the Johnson 
{B, V) fit, with observed colors shifted as described in the text. The lower panel shows the 
(F218W, B) fit. The line along the HB connects ZAHB models of Dorman, Rood, & O'Connell 
(1993), where each actual model is shown as a circle, and effective temperatures of these models 
are shown along the top of the lower panel. The narrow dark lines rising from the ZAHB are 
evolutionary tracks, from the same paper, for HB stars of masses 0.4848, 0.495, 0.510, 0.530, and 
0.540 Mq (from blue to red). The other heavy lines are a zero-age main sequence for stars with 
masses < 1.5Mq, and an isochrone with [Fe/H]= —1.48 and i = 14 Gyr (the age was not fitted to 
the data). 



